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Abstract. A zone diagram is a relatively new concept related to computational 
geometry and fixed point theory. More precisely it is a variation of the notion of the 
Voronoi diagram and formally it is a fixed point of a certain mapping. Neither its 
uniqueness nor its existence are obvious in advance. It has been studied by several 
authors, starting with T. Asano, J. Matousek and T. Tokuyama who considered the 
Euclidean plane with point sites. Computation of zone diagrams is a challenging 
task and in a continuous setting it has been addressed (briefly) only by these authors 
in the above setting. We consider the question of computation of zone diagrams in 
finite dimensional strictly convex spaces with general positively separated closed sites. 
We show that a generalization of the iterative algorithm suggested by these authors 
converges to a double zone diagram, a concept closely related to a zone diagram 
whose existence is known in general. In many cases a zone diagram can be obtained 
from the resulting double zone diagram. We are able to compute approximately 
the corresponding iterations and the resulting (double) zone diagram using a new 
algorithm which enables the computation of Voronoi diagrams in a general setting. 
Unexplained interesting phenomena are discussed too. 

1. Introduction 

1.1. Background. A zone diagram is a relatively new concept related to geometry and 
fixed point theory. In order to understand it better, consider first the more familiar 
concept of the Voronoi diagram. Given a set X, a distance function d, and a collection 
of subsets (Pk)keK in X (called the sites or the generators), we associate with each site 
P k a corresponding Voronoi cell, that is, the set R k of all x E X whose distance to P k 
is not greater than its distance to the union of the other sites Pj, j ^ k. On the other 
hand, in the case of a zone diagram R = (R k )k£K we associate with each site P k the 
set R k of all x G X whose distance to P k is not greater than its distance to the union 
of the other sets Rj, j ^ k. Figures 1 and 2 show the Voronoi and zone diagrams, 
respectively, corresponding to the same ten singleton sites in the Euclidean plane. 

At first sight it seems that the definition of a zone diagram is circular, because the 
definition of each R k depends on R k itself via the definition of the other cells Rj, j ^ k. 
On second thought, we see that, in fact, a zone diagram is defined to be a fixed point 
of a certain mapping (called the Dom mapping), that is, a solution of the fixed point 
equation R = Dom(i2). While the Voronoi diagram is explicitly defined and, hence, its 
existence and uniqueness are obvious, neither the existence nor the uniqueness of a zone 
diagram are obvious in advance. In addition, even if some existence (or uniqueness) 
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Figure 1. The Voronoi diagram Figure 2. The zone diagram of 

of 10 point sites in a square in the the same 10 points as in Figure 1. 

Euclidean plane. 

results are proved, one still faces the problem of finding algorithms for computing zone 
diagrams. 

The concept of a zone diagram was first defined and studied by T. Asano, J. Matousek 
and T. Tokuyama in [2, 3] in the case where (X, d) was the Euclidean plane, each site Pk 
was a single point, and all these (finitely many) points were different. They proved the 
existence and uniqueness of a zone diagram in this case, and also suggested a natural 
iterative algorithm for approximating it. Their proofs and methods rely heavily on the 
above setting. 

Later, more general existence results were obtained. For instance, in [22] it was 
shown that a zone diagram of two general sites in any metric space always exists. 
In fact, the proof holds in a more general setting called m-spaces, in which X is an 
arbitrary nonempty set and the "distance" function d should only satisfy the condition 
d(x,x) < d(x,y) Wx,y 6 X and can take any value in the interval [—00,00]. Simple 
examples given there show that in general uniqueness of the zone diagram does not 
necessarily hold. In [15] it was shown, in particular, that a zone diagram of any finite 
number of sites exists, assuming these sites are compact and positively separated (that 
is, there is a positive lower bound on the distance between any pair of them) and that 
they are located in a large compact subset of lR n with a strictly convex norm (the exact 
result is more general). Another existence (and also uniqueness) result was discussed 
in [12]. Here the setting was the n-dimensional Euclidean space R n , or, more generally, 
finite dimensional normed spaces which are both strictly convex and smooth [12]. The 
sites were again positively separated. 

The existence of a zone diagram is not known in general. However, as proved in [22] , 
the existence of a closely related concept called a double zone diagram is indeed known 
in general (in any m-space, for infinitely many arbitrary sites). Double zone diagrams 
were introduced in [22] as a means for proving the existence of a zone diagram, partly 
because any zone diagram is a double zone diagram. Their existence was shown in a 
nonconstructive way (using the Knaster-Tarski fixed point theorem [14, 25]), and, as a 
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matter of fact, it was shown that there exist a least and a greatest double zone diagrams. 
The importance of double zone diagrams to the computation of zone diagrams will 
become clear later. Zone diagrams are also related to other concepts, such as the 
concepts of trisectors [1, 2, 4, 7], /c-sectors and /c-gradations [10], and territory diagrams 
[8]- 

One of the main challenges regarding zone diagrams is their computation. In a 
continuous setting this task has been addressed so far only by Asano, Matousek, and 
Tokuyama [3] in the case of the Euclidean plane with finite number of distinct points. 
The discussion about an actual method for approximate computing was very brief with 
almost no theoretical or practical details. In spite of the fact that several pictures were 
presented in [3], the only clue regarding the method of computation was mentioned 
briefly at the bottom of page 1184. It was written there that one can use convex 
polygons with many sides for approximating the components of these sequences, and 
that each iteration is computationally demanding. This (approximate) computation 
was restricted to the above setting. As will be explained in Section 3, one of the main 
reasons for this restriction was the lack of a method for computing Voronoi diagrams 
in a general setting. Another reason is the lack of any known representation of the 
boundaries of the involved regions (conjectured to be non-algebraic in many cases). 

In a discrete setting (X is a finite set of points) there has been a limited discussion in 
the case of two sites: a point and a "line segment" or a "curve" in the digital Euclidean 
plane [6] (it was claimed there on the last page that the algorithm can be generalized to 
finitely many line segments but no details were given) and two finite sites in any finite 
m-space [22] . The actual computation in the first case is restricted to the above setting 
and no convergence analysis was given, and there were only clues regarding the time 
complexity. In the second case the actual computation is by brute force. In addition, 
it is not clear in which sense the resulting discrete zone diagrams approximate the 
continuous ones, although intuitively they may approximate them with respect to the 
Hausdorff distance. The computation of double zone diagrams was discussed briefly in 
[22] only in a discrete setting (finitely many finite sites in a finite m-space). 

1.2. Contribution of this paper. We consider the question of computing zone di- 
agrams in finite dimensional strictly convex spaces with general positively separated 
closed sites. We show that a generalization of the iterative algorithm suggested by 
Asano, Matousek, and Tokuyama converges to the least and the greatest double zone 
diagrams. In various cases a zone diagram can be obtained from the resulting double 
zone diagrams, and as a matter of fact, in many of these cases (such as the Euclidean 
case) the limit is the unique zone diagram. We are able to compute approximately the 
corresponding iterations and the resulting (double) zone diagram using a new algorithm 
for computing Voronoi diagrams [19, 20] which enables their approximate computation 
in any normed space (possibly infinite dimensional) and with sites of a general form. 
Many pictures of (double) zone diagrams or their approximating sets produced by the 
method described here are given. It is interesting to note that the main ingredient in 
the computation of each iteration, namely the algorithm for computing Voronoi dia- 
grams, enables the computation of each cell independently of the other ones and hence 
it and the above iterative algorithm can be easily parallelized. The Voronoi algorithm 
[19] can also help regarding the computation of /c-sectors [10] and territory diagrams 
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[8], but these issues will not be dealt with here (see [10, p. 719] for a short related 
discussion). 

1.3. The structure of this paper. In Section 2 we present the basic definitions and 
notation. In Section 3 we discuss qualitatively the iterative algorithm. In Section 4 we 
present the main convergence result. In Section 5 we explain how the corresponding 
iterations are actually computed and we also we briefly discuss the time complexity of 
the algorithm, as well as other topics. The proof of the main result is given in Section 
6. We conclude the paper with Section 7, which contains several interesting questions 
and unexplained phenomena. 

2. Notation and definitions 

In this section we present our notation and basic definitions. We consider a nonempty 
convex subset X in some normed space (X, | • |). We assume that X is not a singleton, 
for otherwise everything is trivial. We note that the definitions given below can be 
generalized almost word by word to the case where (X, d) is a general metric space 
and sometimes even beyond (m-spaces). We denote by A, Int(A) and OA the closure, 
interior, and the boundary, respectively, of the set A . 

Throughout the text we will make use of tuples, the components of which are sets 
(which are subsets of the given set X). Every operation or relation between such 
tuples, or on a single tuple, is done component-wise. Hence, for example, if K ^ 
is a set of indices, and if R = (Rk)k<=K and S = (Sk)keK are two tuples of sets, then 
Rf]S = (Rkf]Sk)keK, R = (Rk)k&K, and R C S means Rk C Sk for each k G K. 
When R is a tuple, the notation (R) k is the fc-th component of R, i.e, (R)k = Rk- 

Definition 2.1. Given two nonempty sets P,ACX, the dominance region dom(P, A) 
of P with respect to A is the set of all x G X the distance of which to P is not greater 
than its distance to A, that is, 

dom(P, A) = {x G X : d(x, P) < d(x, A)}. 

Here d(x, A) = mi{d(x, a) : a G A}. 

Definition 2.2. Let K be a set of at least 2 elements (indices), possibly infinite. Given 
a tuple (Pk)keK of nonempty subsets Pk C X , called the generators or the sites, the 
Voronoi diagram induced by this tuple is the tuple (Rk)keK of nonempty subsets Rk C 
X , such that for all k G K , 

R k = dom(P fe , \JPj) = {xeX: d{x,P k ) < d{x,P 6 ) Vj ^ k, j G K}. 

In other words, each R k , called a Voronoi cell or a Voronoi region, is the set of all 
x G X whose distance to the site P k is not greater than its distance to any other site 
Pj, J ± k. 

Definition 2.3. Let K be a set of at least 2 elements (indices), possibly infinite. Given 
a tuple (Pk)keK of nonempty subsets P k C X , a zone diagram with respect to that tuple 
is a tuple R = (Rk)keK of nonempty subsets Rk C X such that 

R k = dom(P fc , (J Rj) Mk G K. 
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In other words, if we define X k = {C : P k C C C X}, then a zone diagram is a fixed 

point of the mapping Dom : Yl Xk — > Yl X k , defined by 

keK kaK 

Dom(R) = {dom{P k , (J Rj)) k eK. (1) 

a tuple R = (Rk)keK is called a double zone diagram if it is a fixed point of the second 
iteration Dom o Dom, i.e., R = Dom 2 (i?). 

Note that if R is a zone diagram, then it is a double zone diagram as implied by 
iterating the equation R = Dom(i?) with Dom. We now recall the definition of strictly 
convex spaces. 

Definition 2.4. A normed space (X, \ • |) is said to be strictly convex if \x + y\ < 
\x\ + \y\ whenever x and y are any elements in X which are not on the same ray, i.e., 
i^O,?/^ 0, and x/\x\ ^ y/\y\- 

A well known equivalent condition for strict convexity is that the unit sphere of 
the space does not contain line segments. Equivalently, for all x, y G X satisfying 
\x\ — \y\ — 1 and x ^ y, the inequality \(x + y)/2\ < 1 holds. The ci-dimensional 
Euclidean space IR d , or more generally, inner product spaces, and the sequence spaces 
£ p , p G (1, oo), are all examples of strictly convex spaces. On the other hand, M rf with 
the £i or norms are typical examples of spaces which are not strictly convex. See 
[16, 18] for more information regarding strictly convex spaces. 

We finish this section by recalling the definition of the Hausdorff distance. 

Definition 2.5. Given two nonempty sets Ai,A 2 C X, the Hausdorff distance between 
them is defined by 

D(A 1 ,A 2 ) = max{ sup d(a\, A 2 ), sup d(a 2 ,A 1 )}. 

Recall that the Hausdorff distance is different from the usual distance between two 
sets which is defined by d(A ± , A 2 ) = mi{d(ai, a 2 ) : a± G A ± ,a 2 G A 2 }. 

3. A QUALITATIVE DESCRIPTION OF THE ALGORITHM 

As mentioned before, a tuple R = (Rk)keK of cells (subsets) is a zone diagram if 
it satisfies the fixed point equation R = Dom(i?). A common and natural approach 
in fixed point theory for the computation of a fixed point of a given mapping / is by 
iterations [9, 13]. One starts with some point yo in the space Y on which / is defined, 
and starts iterating /. A sequence y\ = f(yo), ■ ■ ■ ,y n +i = /(2/n), • • • is generated, and 
one hopes that it converges in some sense to a fixed point of /. 

Returning to our setting, the given mapping is the Dom mapping. The given space 
on which it is defined is Y = YlkeK -^k, where = {C : P k C C C X}. In other 
words, Y is the collection of all tuples (vectors) whose k-th component is a subset C of 
the given world X such that C contains the site Pk- As a result of the above, a natural 
choice for the starting point yo G Y is the collection (Pk)keK of the given sites. 

It turns out that it is convenient to denote := (Pk)keK = yo, CX°) := Dom(/(°>), 
and define inductively the inner sequence 

j(n) ;= Dom(0 (n " 1) ) and the outer one 

(n) ._ Dom (/W) f or 

any integer n > 1. The reason behind this choice is explained 
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Figure 3. The boundaries of 

/( 2 ) and 0^ in a square in 

tm>2 i, 



The sites are as in Fig- 



ure 2. 




Figure 4. The setting of Figure 
3, but with 1^ and (3) . 



in the following simple lemma, which was essentially observed in [3] and a little bit in 
[20, 22]. 

Lemma 3.1. Let (Pk)keK be a tuple of nonempty and closed subsets in X. 

(a) Dom is antimonotone, i.e., Dom(R) C Dom(S') whenever S C R; Dom 2 is mono- 
tone, that is, RC S => Dom 2 (i?) C Dom 2 (5). 

(b) The sequence (J^)^L is increasing, (O^ n ^)^ is decreasing and we have 

oo oo 

C |J Dom 2 \P k ) keK C f| Dom 2 ^ +1 (P fc )^ C (2) 

7=0 7=0 

for each nonnegative n. In addition, if R is a zone or double zone diagram in X, 
then 

oo oo 

C |J Dom 2 ~<(P k ) keK C P C p| Dom^ 1 ^)^ C > 0. (3) 

7=0 7=0 

Figures 3 and 4 show P 2 -*, O^ 2 ^ and I^ 3 \ respectively, in a square in (M 2 , £2)- The 
sites are as in Figure 2. For another example, see Figure 11. 

From the lemma we see that the goal of I^ 71 ' is to approximate the zone diagram 
from below, and the goal of 0^ n ' is to approximate it from above. However, at this 
stage several difficulties arise. First, it is not clear that these sequences converge, and 
if they do converge, it is not clear that both of them converge to the same limit and 
that the limit of at least one of them is indeed a zone diagram. 

Second, even if a convergence result is obtained, one faces the problem of the actual 
computation of I^ 71 ' and O^ 1 ' for n > 1. Even for n = this is not a simple task since 
= Dom(I^) = (dom(Pfc, LLy/c Pj))k^K is the Voronoi diagram of the sites (Pk)keKi 
and hence, for sites of a general form, or sites which possibly form a degenerated 
configuration, or for space with a general norm, one has to find algorithms which 
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enable the computation of Voronoi diagrams in such a setting, and most of the familiar 
algorithms for computing Voronoi diagrams (see e.g., [5, 17, 24] for some reviews) are 
not able to achieve this task (or, in some cases, they are too complicated or too slow). 
Even if one restricts oneself to the familiar case of the Euclidean plane with point sites, 
one still faces problems starting from the iteration n — 1. This is because one has to 
know how to compute the components of Dom(S') for a given tuple S = (Sk)k^K, and 
hence one has to know a method for computing dom(Pfc, IJ^fc i- e -> the Voronoi 
cell of Pk with respect to the set {Jj^ k Sj. Unfortunately, when S = P n \ n > 1, or 
S = 0^ n \ n > 0, the components of S are general sets, and again, it seems that most 
of the familiar algorithms for computing Voronoi diagrams are not helpful here. 

The sequences (J^ ra ^)^ and (O^)^ =0 were introduced in [3] in the case of the 
Euclidean plane with finitely many distinct point sites. It was shown that IJ^Lo = 
n^Lo and that this tuple is the unique zone diagram. Although this is not a 
pure convergence result (no limits) and although no error estimates for the level of 
approximation were given, one still has the plausible result that (I^)^ =0 increases to 
a zone diagram and (O^ n ^)^L decreases to it. As for the actual computation of 1^ 
and 0^ n \ it was already mentioned in the introduction that very few details were given 
there. 



4. The convergence theorem 

As mentioned in Section 3, it was proved in [3] that in the case of the Euclidean plane 
with point sites one has the equality m := IJ^Lo = C\n=o := and m = M is 
the unique zone diagram. As Theorem 4.1 below shows, in the more general setting of 
strictly convex spaces with possibly infinitely many general positively separated sites 
one has the weaker result that m = Dom(M) and M = Dom(m), and both m and 
M are double zone diagrams rather than zone diagrams (and they are not necessarily 
equal). In addition, Theorem 4.1 also discusses another way to obtain a zone diagram 
in the special case of two sites. 

A corollary to the theorem ensures that when the world X is compact, then (I^)^ =0 
and (O( n ))^ converge to m and M respectively with respect to the Hausdorff distance. 
As a result, since in practice X is taken to be compact, if 1^ and are shown 
experimentally to be almost the same for some n, then one has a good approximation 
to both m and M. Because in this case it is known that a zone diagram exists [15], 
one has a good approximation to this zone diagram, and in fact this also shows that 
probably this zone diagram is unique and coincides with both m and M. In all of the 
pictures presented in this paper this is what we obtained. Another corollary to the 
above ensures that whenever the least and the greatest double zone diagrams coincide, 
then m = M and they both coincide with the unique zone diagram. This is true in 
particular in the Euclidean case. 

Theorem 4.1. Let X be a convex set in a finite dimensional strictly convex space and 
let P = (Pk)keK be a tuple of nonempty closed sets in X with the property that 



mi{d(P k ,P j ):j,keK,j^k}>0. 
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Let 



oo oo 



oo oo 



M = P| Dom 27+1 (P) = p| (n) , 



m — 



(J Dom 27 (P) = (J /("). 



7=0 n=0 



7=0 n=0 



Then 



m = Dom(M) and M = Dom(m), 

andm and M are, respectively, the least and greatest double zone diagrams. If \K\ = 2, 
then by letting m = (mi,m 2 ) and M = (Mi,M 2 ) 7 each of the pairs (mi,M 2 ) and 
(Mi,m 2 ) is a zone diagram. 

Corollary 4.2. Under the setting of Theorem if X is compact, then by letting 
M = (Mfc)fc e x and m = (mk)keK we have 



with respect to the Hausdorff distance. 

Corollary 4.3. Under the setting of Theorem 4-1, if it is known that the least and the 
greatest double zone diagrams coincide, then m = M and they both coincide with the 
unique zone diagram (even if X is not compact). In particular this is true when the 
normed space is Euclidean. 

Proof. By [22, Corollary 6.2] (whose proof is elementary) if the least and the greatest 
double zone diagram coincide, then, without any restriction on the sites or the space, 
there exists a unique zone diagram which is equal to both of them. Since in the proof 
of [12, Theorem 1.1] it was proved that the least double zone diagram coincides with 
the greatest one when the normed space is Euclidean, the assertion follows. □ 

Example 4.4. An illustration of Theorem 4.1 for the case where m ^ M is given in 
Figures 6-8. Here the sites are Pi = {(0,1)}, P 2 = {(0,-1)} in a rectangle in M 2 , 
but the norm is ||(xi, x 2 )|| = 5^J a 2 \xi\ 2 + |a; 2 | 2 + (1 — + (1 — 5)|x 2 |, where 

a = 5 = 0.1. This norm is strictly convex but not smooth (the unit sphere contains 
points which have several supporting lines). A closely related example was discussed in 
[12, Section 5] in the context of non-uniqueness of zone diagrams (double zone diagrams 
were not mentioned) and it actually follows from the discussion presented there that 
indeed m ^ M in this case. Two different zone diagrams can be obtained now, by 
taking either the first component of M together with the second of m or the second 
component of M together with the first of m. 

For producing the pictures (see also Section 5) the endpoints of 4000 rays emanating 
from each site were computed in the Voronoi algorithm stage. However, it can be seen 
that still some parts of the boundaries of M are not full. The reason is that many rays 
should be produced in a very small angle (i.e., a very small part of the unit sphere) 
because of the location of the site with respect to these boundary parts. A simple way 
to overcome the problem (after detected) is to distribute the rays in the unit sphere in 
a non-uniform way such that many rays will emanate in the corresponding angle but 
relatively few in other angles, in contrast with the current way of producing the rays 
in a uniform way. One can also connect neighbor endpoints by line segments. 

The proof of Theorem 4.1 is quite long and technical and it can be found in Section 6 
together with the proof of Corollary 4.2 which follows from a quite general argument not 



M k = lim (Dom 27+1 (P)) fc and m k = lim (Dom 27 (P)) fc \/k e K 




(4) 
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Figure 5. The components of 
the greatest double zone diagram 
M from Example 4.4. 




Figure 7. J (2) and (2) for the 
double zone diagrams from Ex- 
ample 4.4. 



Figure 6. The components of 
the least double zone diagram m 
from Example 4.4. 




Figure 8. J (3) and (3) for the 
double zone diagrams from Ex- 
ample 4.4. 



actually related to zone diagrams. Some of the ideas used in the proof of Theorem 4.1 
are partly inspired from the ones given in [3, Lemma 3.1, Lemma 5.1]. However, the 
generalization of the arguments given there to our setting is definitely not immediate 
and one has to pay attention to certain subtle points, among them the verification of 
the equality 

oo oo 

jj f|(Dom 2 ^ +1 (P fc ) fe6i ,) j = p| UtDom^PfcW),- 

jjLk 7=1 7=1 j^k 

which is based on some non-immediate properties of dom and Dom. 

It is interesting to note that Theorem 4.1 is also related to [10, Proposition 5] in the 
sense that in both cases one considers an iterative algorithm and the setting is closed 
sites in a finite dimensional strictly convex space. The proofs of these two claims 
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have some similarities: for instance, [10, Lemma 6] is essentially Proposition 6.5 in 
Section 6. However, the proofs of these two claims are different in the strategy and 
the level of detailing (and historically they have been established independently). In 
general the proofs of [10, Proposition 5] and Theorem 4.1 here are different since the 
involved mappings are not the same and hence one needs to find a strategy which is 
appropriate for each case separately. In fact, an important feature used in the proof 
of Theorem 4.1 is the fact that the sites are positively separated (this also allows us 
to consider infinitely many sites) while in [10, Proposition 5] the sites satisfy the more 
general condition of being merely disjoint, but now one must consider only two sites 
(by definition) and the mapping must have finitely many components. 

5. The actual computation of J (n) and (n) , time complexity, and 

CLARIFICATIONS 

In this section we provide details about practical and theoretical aspects of the 
method which allows the approximate computation of the corresponding regions. The 
section has three part: in the first part practical details regarding the method are 
given. In the second part the issue of time complexity is considered. In the last part 
several related issues are clarified. 

5.1. The actual computation of 1^ and O^: As mentioned in Section 3, in order 
to compute the iterations 1^ and one has to know how to compute or approximate 
dominance regions dom(P, A) of general sets P, A. The way we choose to overcome 
this difficulty is to use the approximation algorithm for computing Voronoi diagrams 
of general sites in general normed spaces which was introduced in [19]. In a nutshell, 
for approximating dom(P, A) one uses the fact that dom(P, A) can be represented as 
a union of line segments emanating from the points of P (see Figures 9-10 for an 
illustration). 

First, the world X is assumed to be a large compact subset, e.g., a square or a 
hypercube. Now one approximates P using a finite collection of points; this is always 
possible if P is compact. After that a finite collection of directions (unit vectors) is 
chosen, and one approximates the endpoints of the line segments emanating from the 
points of P in these directions. These endpoints can be approximated to any required 
precision, and then they are stored. At this point dom(P, A) is represented by this 
collection of endpoints and it is regarded as computed. 

For computing 7^ and one computes their corresponding components iter- 
atively: the components are dom(P fc , A k ) where k runs over all the (finitely many 
in practice) indices in K and A k depends on the iteration and on k. In general, 
A k = {j^ k {I {n) ) k or A k = U^fc(° (n) )fc where ( /(n) )fc is the k - th component of I (n \ 
For instance, the components of are dom(P k , A k ) where A k = {_}j^ k {0^) k and 

(O^)k = dom(Pfc, (Jj^jfc Pj) i s known, as a collection of endpoints, from a previous 
computation. An illustration of the process is given in Figures 3,4,11, and 12. 

An important parameter which determines the level of approximation is the number 
of line segments (or, actually, the density of the corresponding unit vectors in the unit 
sphere of the space) used in the construction of the dominance regions. In Figures 3 
and 4 about 1600 line segments emanate from each site, and only the endpoints are 
displayed. In Figures 8-fig:ZD-Iteration2-StrangeNorm-2Sites-0002 4000 rays emanate 
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Figure 9. The Voronoi dia- 
gram of 4 sites in (M. 2 ,£ p ), p = 
2.718281828, each consists of 3 
points; 80 line segments emanate 
from each point of each site. 




Figure 10. The zone diagram 
(and hence the double zone dia- 
gram) of the sites of Figure 9. 



from each site and only endpoints are displayed. In Figures 9-10 about 80 line segments 
are used, and the whole line segment is displayed. 

5.2. Time complexity: As for the time complexity of the algorithm, it is determined 
by the time complexity of dom(Pfc, A^) for each k which is a particular case of the 
algorithm for computing Voronoi diagrams described in [19, 20]. Starting from iteration 
n — 1, the number of endpoints in is CiC 2 ~ 1 G(\K\ — 1) where C\ is a universal 
positive constant depending on the norm, C 2 is a positive constant depending on the 
level of approximation of the unit sphere of the space (roughly speaking, it is a bound 
on the distance between a point on the unit sphere and some unit vector from the 
list of unit vectors created by the user), d is the dimension, G is a bound on the 
number of points in (for each k), C 1 C , 2 _1 G is the number of line segments emanating 
from the points of and \K\ is the number of sites. For determining an endpoint 
up to some error parameter one makes C3 1^4^ | comparisons where C3 is a positive 
constant depending on the level of approximation. Hence the time complexity for the 
entire number of components is OiC^CxC^^G^^ 2 ). Note that here G and \K\ are 
parameters depending on the input, G\ and d are global parameters, and C2 and C3 
are parameters depending on the user (the level of approximation of the computation). 
Under some assumptions, e.g., that the sites are points which are uniformly distributed, 
it is possible to reduce significantly the number of calculations. 

5.3. Clarifications: We finish this section with a discussion on certain issues which 
perhaps have not been so clear so far. It should be emphasized that the goal of 
this section is to describe schematically a practical way for approximating the regions 
which appear in each iteration and to evaluate roughly the number of calculations. The 
fact that the algorithm can approximate a given dominance region up to any desired 
precision is a simple consequence of the stability of the algorithm mentioned in [19]. 
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However, the proof of this stability property is not obvious and it will be discussed 
elsewhere in a paper devoted to this algorithm. Its proof is in the spirit of the proof of 
the geometric stability of Voronoi diagrams with respect to small changes of the sites 
[21]. 

Roughly speaking, if the users want to approximate (in terms of the Hausdorff dis- 
tance) a given region up to some error parameter e, then they need to approximate well 
enough the sites, to choose in advance enough approximating rays (this is determined 
by the error parameter related to the unit sphere), and to fix a small enough error 
parameter for the endpoints of the rays. If the users decide in advance how many iter- 
ations they want to perform for approximating the (double) zone diagram and what is 
the level of approximation of the regions in the final iteration, then they can estimate 
in advance the number of calculations by iterative "reverse engineering" : using the tar- 
get error parameter e, one estimates the error parameters needed as input for the final 
iteration, and from them the error parameters needed for the previous iteration, and 
so on, until the initial iteration. This gives an estimate on the initial error parameters. 

In the above description the number of iterations n was chosen by the user but there 
was no guarantee that the real (double) zone diagram will be approximated well by 
j(«) Qr (n)_ However, since the algorithm converges to the double zone diagram by 
Theorem 4.1, then given e > 0, there is a number no, depending only on e, such that 
for any n > no, the regions of any iteration 1^ and will be at Hausdorff distance 
of at most e from the limit regions. Hence, if one can compute no(e), then one can 
know in advance the corresponding needed initial error parameters. Unfortunately, it 
is not clear how to estimate no(e) and this is a major open problem. 



6. Proof of the main result 

In this section we present the proof of Theorem 4.1 and the corollaries following it. 
We start with several lemmas describing simple properties of dom and Dom. 

Lemma 6.1. Let (X,d) be a metric space. Then 

(a) d(x, U 7 er ^-7) = inf{d(x, Ay) : 7 G T} for any x G X and any collection {Ay} 7£ r 
of nonempty subsets in X . 

(b) dom(P, IJ 7 er ^7) = Pl 7 er dom(P, Ay) for any collection {Ay} 7e r of nonempty sub- 
sets in X and any PCX nonempty. 

(c) dom([J 7gr P 7 , A) C |J 7gr dom(P 7 , A) for any collection {P 7 } 7e r of nonempty sub- 
sets in X and any A C X nonempty. If, in addition, for each x G X the distance 
d(x,\J ieF P 7 ) is attained at some P 7o; i.e., d(x, |J 7e r ^7) = d(x,P 10 ) for some 
70 G T ; then the equality dom(|J 7gr P 7 , A) = IJ 7 er dom(P 7 , A) holds. 

(d) Dom(lJ 7gr P 7 ) = H 7 er Dom(P 7 ) for any collection {P 7 } 7er of tuples, where each 
K 1 is indexed by the same set of indices K, i.e., K 7 = (Rj)keK- Here Dom is 
defined with respect to some tuple (Pk)keK of nonempty subsets in X . 

Proof, (a) Let x G X, a := d(x, |J 7gr Ay) and (3 := inf{<i(x, Ay) : 7 G T}. Then 
a < d(x, Ay) for all 7 G T by the definition of a, so a < (5. If a < (3, then there 
is y G U 7 er ^7 sucn ^ na ^ ^( x > v) < Pi anc ^ smce V £ Ay for some 7 G T, we have 
d(x, Ay) < d(x, y) < (3, a contradiction with the definition of (3. 
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(b) By the antimonotonicity of dom(P, •) we have dom(P, lj igr A) C dom(P, Ay) for 
any 7 G Y. Consequently, we have dom(P, (J 7 er ^7) — fl 7 gr dom(P, Ay). Con- 
versely, suppose that x G f] 7 er dom(P, Ay). If x ^ dom(P, (J 76r Ay), then there 
is y G U 76r Ay such that d(x,P) > d(x,y). But y G Ay for some 7, and 
x G dom(P, Ay), so d(x,P) < d(x, Ay) < d(x,y), a contradiction. 

(c) By monotonicity of dom(-,A), i.e., dom(P, A) < dom(Q, A) whenever P C Q 
(which is easily verified), we have dom(P 7 , A) C dom(|Jj gr Pj, A) for any 7. As a re- 
sult, we have IJ 7 er dom(P 7 , A) C dom(|J 7Gr P 7 , A). Now suppose that d(x, IJ 7e r ^7) = 
d(x,P l0 ) for some 70 G T. Therefore, if x G dom(lJ 7gr P 7 , A), then d(x,P l0 ) = 
d ( x ^er P i) ^ d(x,A), i.e., x G dom(P 7o ,A) C \J yer dom(P 7 , A). 

(d) Since the union and intersection are taken component-wise, this follows from the 
definition of the Dom mapping (equation (1)) and part (b): 

Dom(|J W) = (dom(P fc , \J(\J R]))) keK = (dom(P fc , \J({J R]))) keK 
7 er j^k 7 er 7 er 

= (f| dom(P fe , |J R])) keK = f) Dom(P^). 

7 er j^k 7 er 

□ 

Lemma 6.2. Let (P k ) k( zK be a tuple of nonempty and closed subsets in X. 

(a) Dom is antimonotone, i.e., Dom(P) C Dom(S') whenever S C R; Dom 2 is mono- 
tone, that is, RCS Dom 2 (P) C Dom 2 (S). 

(b) The sequence (I^)™ =0 is increasing, (O^)^ =0 is decreasing and we have 

00 00 
C |J Dom 2 ^(P fc )^ C f| Dom 2 ^ +1 (P k ) keK C (5) 

7=0 7=0 

for each nonnegative n. In addition, if R is a zone or double zone diagram in X, 
then 

00 00 
C |J Dom 2 ^(P fe ) fceX C R C f| Dom 2 ^\P k ) keK C Vn > 0. (6) 

7=0 7=0 

Proof, (a) The antimonotonicity of Dom is a simple consequence of the easily verified 
fact that dom(P, B) C dom(P, A) for any nonempty subset P and any nonempty 
subsets A, B satisfying A C B. The composition of two antimonotone mapping is 
monotone, and hence Dom 2 is monotone. 

(b) Let {X) k& x be the tuple whose components are X. An immediate check shows 
that P k = dom(P fc ,X) and Dom(X) k< z K = (P k ) ke K- From the inclusion (P k ) k ^K Q 
(X) keK , the definition of the space on which Dom acts, and the monotonicity of 
Dom 2 we have 

(Pk)keK C Dom 2 (P k ) keK C Dom 2 (X) keK = Dom(P k ) keK C (X) keK . 

As a result, by iterating with Dom 2 , using its monotonicity, and using the definition 
of and we see that (I^)'^ ) =0 is increasing and (O^)^L is decreasing. 
Hence is contained in {J e=Q I^ and f] i=0 O^ is contained in for each n. 
Finally, after iterating the inclusion (P k ) k& K Q R Q (X) ke K using Dom 2 we obtain 
(6) when R is a double zone diagram. But a zone diagram is a double zone diagram 
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[Dom(Dom(P)) = Dom(P) = R whenever R = Dom(P)], so (6) is true in this case 
too. 

□ 

Lemma 6.3. Let (X, d) be a metric space and let P,ACX be nonempty. Then 

(a) dom(P, A) is a closed set. 

(b) dom(P, A) = dom(P, A) = dom(P, A) = dom(P, ~A\ 

(c) For a tuple R = (R k ) k eK of nonempty subsets, let R = (R k ) ke K- Then Dom(P) = 
Dom(P) = Dom(P) (with respect to a given tuple (P k ) k ei<)- 

Proof, (a) If x = lim^oo x n and x n G dom(P, A), then d(x n ,P) < d(x n ,A), and this 
inequality is preserved in the limit because the function x h- > d(x, P) — d(x, A) is 
continuous (with respect to the topology induced on X by d). 

(b) This follows from the general fact that d(x, B) = d(x, B) for each x G X and each 
PCX. 

(c) By part (b) and by Lemma 6.1(b), we have 

dom(P fe , |J R~) = p| dom(P fe , R~) = f] dom(P fe , R 3 ) = dom(P fc , |J Rj). 

jy^k jy^k jyik j^k 

Thus Dom(P) = {dom{P k ,{J^ k R~)) k&K = (dom(P fe , (J^ Rj))keK = Dom(P). 
Finally by part (a) and by definition, we have 

Dom(P) = (dom(P fc , \J Pj)) keK = (dom(P fe , |J Pj)) keK = Dom(P). 

j^k j+k 

□ 

Proposition 6.4. Let X be a convex subset of a finite dimensional strictly convex 
normed space (X, | • |). Let A C X be nonempty and closed. Let p,z e X and suppose 
that d(z, p) < d(z, A) and p ^ A. Then d(x, p) < d(x, A) for all x e [p, z) . 

Proof. We can assume that z ^ p, because otherwise the assertion is obvious (void). 
Fix x G [p, z) and let a x G A be chosen such that d(x, A) = d(x, a x ). There exists such 
an a x by compactness since the space is finite dimensional and A is closed. Since p ^ A 
and d(z,p) < d(z,A) it follows that z,p and a x are all different. 

We distinguish between several cases. Suppose first that p, z, a x are on the same 
line. If p is between z and a x , then obviously d(x,p) < d(x,a x ) = d(x,A). The case 
where a x is between z and p is impossible, since then d(z,a x ) < d(z,p), contradicting 
d(z,p) < d(z,A). If z is between p and a x , then d(z,a x ) < d(z,x) + d(x,a x ). This 
inequality is true also if p, z and a x are not on the same line, by the strict convexity of 
the norm. Hence, in these two latter cases we have 

d(z, x) + d(x, p) = d(z, p) < d(z, A) < d(z, a x ) < d(z, x) + d(x, a x ), 

so d(x,p) < d(x, a x ) = d(x, A). □ 

Proposition 6.5. Let X be a convex subset of a finite dimensional strictly convex 
normed space (X,\ ■ |). Let P,ACXbe nonempty and closed, and suppose that 
PC\A = ®. Then 

dom(P,A) = {x G X : d{x,P) < d{x,A)}. 
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Proof. From the inclusion {i6l: d(x,P) < d(x,A)} C dom(P, A) and the fact that 
dom(P, A) is closed (Lemma 6.3(a)) we have the inclusion {x G X : d(x, P) < d(x, A)} C 
dom(P, A). For the second inclusion, let z G dom(P, A). If d(z,P) < d(z,A), then 
obviously the point z is in the set {x & X : d(x,P) < d(x,A)}. Now suppose that 
d(z ) P) = d(z,A). By our assumption there is p G P such that d(z,P) = d(z,p), 
and since p G P, it follows that p ^ A. Since A is closed we have d(p, A) > 0. 
Hence if z = p, then z G {re G X : d(x, P) < d(x, A)}. If z ^ p, then [p, z) 7^ 0, and 
Proposition 6.4 implies that any x G [p, z) (arbitrary close to z) satisfies the inequality 
d(x,P) < d(x,p) < d(x,A). Thus again z G {x G X : d(x,P) < d(x, A)}. 

□ 

The following lemmas are somewhat independent of the main concepts of the paper 
but they are needed later. 

Lemma 6.6. Let X be a closed subset in a normed space (X, \ ■ |) and let (A J )^ =1 be 
a decreasing sequence of subsets of X such that A : = H^Li ^7 ^ s nonempty. If A 1 is 

closed for each 7 and (X, | • |) is finite dimensional or X is compact, then d(x,A) = 
lim^oo d(x, Ay) for each x G X . 

Proof. Let x G X and assume that (X,\ ■ |) is finite dimensional or X is compact. In 
this case compactness arguments yield that the distance between x and any nonempty 
closed subset is attained. Therefore for each 7 there exists x 7 G Ay such that d(x, Ay) = 
d(x,Xj). Since (A y )^ =1 is decreasing and Ay is closed for each 7, it follows that 
any limit point of the sequence (:c 7 ) 76 n is in A 1 for each 7 and hence in A. Since 
(A y ) 7 < L 1 is decreasing, the sequence (d(x, x 7 )) 76 pj is increasing and its limit is bounded 
by d(:r,A). Hence (a; 7 ) 7g N is bounded and there exists Xoo G A and a subsequence 
such that a;^ = lim^oo x 7 . From the continuity of the distance function we deduce 
that d(x,x 00 ) = lim^oo d(x, x lp ) < d(x,A). But d(x,A) < d^x^x^) by definition, so 
lim^oo d(x, Ay) = lim^oo d(x, A lfj ) = d(x, Xoo) = d(x, A). □ 

Lemma 6.7. Let (X,d) be a compact metric space. Let (A n )^ =1 be a decreasing se- 
quence of nonempty closed sets of X , and let (P n )$£Li be an increasing sequence of 
nonempty sets of X . Let A = H^Li A n an d B = U^Li B n - Then A = lim^oo A n and 
B = lim^oo B n with respect to the Hausdorff distance. 

Proof. Suppose by way of contradiction that it is not true that lim^oo D(A n ,A) = 0, 
where D is the Hausdorff distance. Then there exists some e > and a subsequence 
n m G N such that e < D(A nm , A) for each m G N. Since A C A n for each n and because 
X is compact, the definition of the Hausdorff distance implies that for each m there 
exists x m G A nm such that e < d(x m ,A). Let (x mi )i e ^ be a convergent subsequence of 
the sequence (x m )™ =1 , and let x be its limit. Since the sequence (A n ) n is decreasing, it 
follows that x m , G A„ whenever I < V . Therefore x G A, because A„ is closed. 
Thus x G f]^i 4m, • But Pli^i ^n mi = A because the intersection is decreasing. Hence 
= d(x,A) = lim/^oo d(x m ,A) > e, a contradiction. 

Now suppose by way of contradiction that it is not true that lim^oo D(B n , B) = 0. 
Then there exists some e > and a subsequence n m G N such that e < D(B nm , B) for 
each m G N. Since B n C B for each n and X is compact, the definition of the Hausdorff 
distance implies that for each m there exists x rn G B such that e < d(x m , B nm ). Let 
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(xmJjeN be a convergent subsequence of the sequence (x m )™ =1 , and let x be its limit. 
Then d(x,x m ) < e/4 for all / large enough. Since B is closed it follows that x G B, 
so by the definition of B there exists n G N such that d(x,B no ) < e/2. But (B„)^ =1 
is increasing, so d(x,B n ) < e/2 for all n G N such that n > n , and, in particular, 
d(x ) B rimi ) < e/2 for / large enough. Thus 

e < d(x mi ,B nmi ) < d(x mi ,x) + d(x,B nmi ) < 3e/4, 

a contradiction. □ 

In the next claims more subtle properties of dom and Dom will be derived. Lemma 
6.10 and Proposition 6.11 generalize [3, Lemma 3.1] and [3, Lemma 5.1] respectively, 
and are partly inspired by them. 

Lemma 6.8. Let (X,d) be a metric space and let P = (P k ) k <=K be a tuple of nonempty 
subsets of X such that for all k G K , 

r k :=mi{d(P k ,P 1 ):j^k}>0. 

For QCI nonempty and r > 0, let B(Q,r) = {x G X : d(x,Q) < r}. We denote by 
(DomP)fc the k-th component o/DomP. Then: 

(a) B{P k ,r k /A) C (Dom 2 P)fc C (Dom 7 P)fc for any integer 7 > 1 and any k G K. 

(b) If X is a convex subset of a normed space, then d(x, B(Q, r)) + r < d(x, Q) for any 
Q C X nonempty, r > and x B(Q,r). 

(c) If d(x, B(P k ,r)) < d(x,P k ) for any k G K, r > 0, and x P k , then for each 
integer 7 > 2 the components of Dom 7 P are disjoint. Moreover, if X is a convex 
subset of a normed space, then for any j, k G K ,j 7^ k and any 7 > 2 we have 
r k /8 + rj/% < rf((Dom 7 P) fe , (Dom 7 P) i ). 

(d) Suppose that X is a convex subset of a normed space. Given j, k G K , j 7^ k, and 
an integer^ >2, if x G (Dom 7 P)j ; then r k /A< d(x,P k ). 

Proof, (a) Let k G K and suppose that x G B(P k ,r k /A), i.e., d(x,P k ) < r k /4. By 
definition we have (Dom 2 P) fe = dom(P fc , (Jj^fc dom(Pj, (J^ . Pi)), so in order to 
prove that x G (Dom 2 P) fc it suffices to prove that r k /A < d(x,y) for all y G 
U^don^Pj, [Ji^j Pi)- Given y in this union, it is in dom(P,, \J izjtj Pi) for some 
j 7^ k. Hence 

d(y, Pj) < d(y, |J Pi) < d(y, P k ) < d{y, x) + d(x, P k ). 

Therefore, by the definition of r k , 

r k < d{Pj, P k ) < d(Pj, y) + d{y, x) + d(x, P k ) < 2d(x, y) + 2r k /4, 

i.e., d(x, P k ) < r k /A < d(x,y) and the assertion follows. Finally, Lemma 6.2 implies 
that Dom 2 P C Dom 7 P for any integer 7 > 1, and hence (Dom 2 P)fc C (Dom 7 P)fc 
for each k G K. 

(b) Let p G Q. Since x B(Q,r), a simple argument shows that the intersection of 
the compact segment [x,p] with the closed set dB(Q,r) is a nonempty compact 
set, and hence the distance between x and this intersection is attained at some 
point y G [x,p\. It must be that s := r — d(y,p) < 0. If this is not true, then 
for any z in the open ball B(y,s) we have d(z,Q) < d(z,p) < d(z,y) + d{y,p) < 
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s + d(y,p) = r, and hence y is in the interior of B(Q,r), a contradiction. Thus 

d(x, p) = d(x, y) + d(y, p) > d(x, y) + r, so 

d(x, B(Q, r)) = d(x, B(Q,r)) < d(x, y) < d(x, p) - r, 

for any p G Q. As a result, d(x, B(Q,r)) + r < d(x, Q). 

(c) Suppose that we know that the components of Dom 3 P are disjoint. Then for each 
integer 7 > 2 the components of Dom 7 P are disjoint, because Dom 7 P C Dom 3 P 
by Lemma 6.2. We now give a proof for the case 7 = 3. Let k\ 7^ ft2 be two indices 
in K, and assume to the contrary that x G (Dom 3 P)fc 1 f](Dom 3 P) k2 . By definition 
and by part (a) we have 

d(x,P kl ) < d(x, |J (Dom 2 P),) < d(x, (Dom 2 P) fe2 ) 

<d(x,B(P k2 ,r k2 /4))<d(x,P k2 ). 

Hence x ^ P k2 , because otherwise also x G P kl , a contradiction with the assumption 
that < r kl < d(P kl , P k2 ). Thus d(x, B(P k2 ,r k2 /A)) < d(x,P k2 ) by assumption, so 
actually d(x,P kl ) < d(x,P k2 ). In the same way d(x,P k , 2 ) < d(x,P kl ), a contradic- 
tion. 

Finally, suppose that X is a convex subset of a normed space. Let j,k G K be 
different. Since (Dom 7 P)j C (Dom 3 P)j for each integer 7 > 2 and each i G K, it 
suffices to show that r k /8+r-j/8 < d(x, y) for any x G (Dom 3 P)fc and y G (Dom 3 P)j. 
By definition, the triangle inequality and parts (a),(b), 

d(x,P k ) < d(x, y(Dom 2 P)i) < d(x, (Dom 2 P)j) < d(x,B(P j ,r j /A)) 

< d(x, Pj) - r,/4 < d(x, y) + d(y, P 3 ) - r,/4. 

In the same way d(y, Pj) < d(y, x) + d(x, P k ) — r k /A. By adding these two inequal- 
ities we obtain the desired conclusion. 

(d) If x G (Dom 7 (P)) i , then x <£ B(P k ,r k /4) C (Dom 7 (P)) fe by parts (a) and (c). 
Hence r k /A < d(x,P k ) by part (b). 

□ 

Remark 6.9. Lemma 6.8 parts (b),(c),(d) hold with exactly the same proof in any 
geodesic metric space. Recall that a metric space (X, d) is called geodesic if between 
any two points x and y in it there exists an isometry 7 (i.e., 7 preserves distances) from 
the real line segment [0, d(y, x)] to X such that 7(0) = x and 7(<i(?/, x)) = y. Simple 
and familiar examples of geodesic metric spaces are: the Euclidean plane, any convex 
subset of a normed space, spheres, complete Riemannian manifolds [11, pp. 25-28], 
and hyperbolic spaces [23, pp. 538-9]. 

Lemma 6.10. Let (X,d) be a metric space, let P C X be nonempty and suppose that 
{A i \ ( *L l is a family of subsets of X such that H^i A 0- U 



00 00 
dom(P, p| A) = {x G X : d(x, P) < d(x, f| A t )}, 
i=i i=i 



(7) 
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and 

oo 

d(y, p| Ai) = lim sup d(y, A t ) , Vy G X, (8) 



i=i 



dom(P, P| At) = |J dom(P, A t ) . (9) 

i=l i=l 

Proof. By antimonotonicity, dom(P, f) ( *L 1 Aj) D dom(P, A*) for all i, and since the left 

hand side is closed, the inclusion dom(P, H^i A) ^ Ui=i dom(P, Aj) follows. For the 
reverse inclusion, let e > be given and suppose that x G dom(P, Ai). We should 
prove that there are % G N and y G dom(P, Aj) such that x) < e. 

By (7) there is y G X such that d(x,y) < e and r := d(y,f]°l 1 A i ) — d(y,P) > 0, 
and by (8) there is % large enough such that \d(y, Hjli A?) — ^(f> < r /2- Hence 

oo 

%, P) + r/2 < %, p| Aj) < d(y, A t ) + r/2. 
i=i 

Thus d(y, P) < d(y, Ai), and so y G dom(P, A { ). □ 

Proposition 6.11. Let X be a be a convex set in a finite dimensional strictly convex 
space (X, | ■ |) and let (P k ) ke K be a tuple of nonempty closed sets in X such that 
r := inf {d(P k , Pj) : j, k G K, j ^ k} > 0. Then 



Dom(p Dom^' +l (P k ) keK ) = [j Dom^(P k ) keK . 

7=0 7=0 

Proof. For each nonnegative integer 7, let P 7 = Dom 27+1 (P k ) ke K- We claim that 

00 00 

ui>j=niK- ( io ) 

7=1 7=1 jyfe 

Indeed, if x G U^fc D^Li ^J> tnen x e D^Li -^J f° r some i 7^ ^. Since P 7 C U^ fe PJ, 
it follows that P|^=i R ] ^ fl^li ^7> so x G fl^Li #7- 0n tlie other hand, 
let x G D^zi Uj^it R ]- Given 7 G N, Lemma 6.8(c) (note that 27 + 1 > 3 and 
P 7 = Dom 27+1 (P fc ) fce x) implies that P 7 f\R] = for j 7^ i, so there exists exactly one 
index j 7^ k such that x G P 7 . It must be that all these indices coincide. If this is 

not true, then x G R~jf)Rj, for some 7' > 7 and corresponding indices j' 7^ j. But 
R]> Q R]/ by Lemma 6.2, so x G R~jf)R^, a contradiction to Lemma 6.8(c). As a 
result x G H^=i -^J f° r s °me j 7^ fc, and this establishes (10). 

Now fix k G K and let A 7 = {Jj^ k R'J and A = {^ =1 ^7- By Lemma 6.2 we 
have Pj C f^^P 7 , so (10) implies that A ^ |. By Lemma 6.8(d) it follows that 
0<r/A<d{P k ,A). 

By Lemma 6.2 the intersection which defines A is decreasing. By Lemma 6.3(a) and 
Lemma 6.8(c) each A 7 is a closed set, because it is a union of closed and disjoint sets 
with a positive distance r/4 between any two different members in the union. Therefore 
Lemma 6.6 implies that d(x, A) = linx^^ d(x, A 7 ) and hence (8) holds. In addition, 
from the above we see that A is closed and hence we conclude from Proposition 6.5 
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that (7) holds for P = P k . Consequently, we conclude from Lemma 6.2, the definition 
of Dom, (10), and Lemma 6.10 that 



oo 

7 



Dom ( P| Dom 27+1 (P fe ) fce ^ = Dom If] W ) = Dom ( f] R 

VY=0 / \7=0 



= dom [P k ,{Jf]R])} = dom (P k ,f]{jK 



keK ^ ^ 7=M^K / / keK 




|J dom(P fc , (J Pj) = |J Dom(PT) = |J Dom 2 \P k ) keK . 

V7=l j^fc / fe e ^ 7=1 7=0 

□ 

We are now able to prove the main result. 

Proof of Theorem 4.1. Part of the proof is completely general. By Lemma 6.1(d) 
and Lemma 6.3(c) we obtain 



Dom(m) = Dom ( (J Dom 27 (P) j = f] Dom 27+1 (P) = M. (11) 

\7=0 / 7=0 

On the other hand, Proposition 6.11 implies that 



Dom(M) = Dom f] Dom 27+1 (P) J = (J Dom 27 (P) = m. (12) 

\7=0 / 7=0 

Thus m = Dom(M) = Dom(Dom(m)) and M = Dom 2 (M), so both m and M are 
double zone diagrams. From Lemma 6.2 and Lemma 6.3(c) we see that m and M 
are respectively the least and greatest double zone diagrams. Finally, if \K\ = 2, 
then by (11), (12), and the definition of Dom we have (Mi, M2) = Dom(m) = 
((dom(Pi, m 2 ), dom(P 2 , mi)) and (m 1? m 2 ) = (dom(Pi, M 2 ), dom(P 2 , Mi)). Hence mi = 
dom(Pi,M 2 ) and M 2 = dom(P 2 ,mi). But the definition of Dom implies the equality 
Dom(mi,M 2 ) = (dom(Pi, M 2 ), dom(P 2 , mi)). Therefore Dom(mi, M 2 ) = (m 1; M 2 ), 
i.e., (mi,M 2 ) is a zone diagram. By the same way (M 1; m 2 ) is a zone diagram. 

□ 

Proof of Corollary 4.2. By Lemma 6.2 the sequence A y = (Dom 27+1 (P)) fc is de- 
creasing and the sequence B y = (Dom 27 (P)) fc is increasing for each keK. As a 
result, when X is compact (4) follows from Lemma 6.7. □ 

7. Open problems 

We conclude this paper with several interesting open problems and unexplained 
phenomena. First, perhaps the most interesting open problem is to establish error 
estimates for the convergence speed of (I^)^ =0 and (O^)^ =0 . One of the reasons that 
this problem is interesting, is the observed fast convergence of these sequences (usually 
4 iterations suffice). We feel that the approach and error estimates developed in [21] 
may help here. 
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in (IR 2 ,^); the bold line is 
/( 3 ) « « « OW. See 
also [22, Example 2.4]. 



Second, it is interesting to find out whether Theorem 4.1 can be generalized to 
other settings, for instance to all normed spaces, or at least to finite dimensional ones. 
Recently we have made some progress regarding this issue in the particular case of 
infinite dimensional uniformly convex spaces. This issue will be discussed in a future 
version. 

Third, although in general a zone diagram is not necessarily unique, even in the case 
of M 2 with two point sites (see e.g., [12, Sections 1,5] or Example 4.4), and although 
our results are limited to strictly convex spaces, experiments show that the algorithm 
actually converges to a unique zone diagram most of the times also in the case of 
non-strictly convex spaces; see Figures 11-12 for an illustration. 
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